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ABSTRACT 

Matter accreting onto the magnetic poles of a neutron star spreads under gravity towards the 
magnetic equator, burying the polar magnetic field and compressing it into a narrow equatorial belt. 
Steady-state, Grad-Shafranov calculations with a self-consistent mass-flux distribution (and a semi- 
quantitative treatment of Ohmic diffusion) show that, for > 10~^ Mq, the maximum field strength 
and latitudinal half- width of the equatorial magnetic belt are fima? 



and Ae* = max[3°(Ma/10-''MQ)-i-5,3°(Ma/10-''Mo)0-5(Ma/10 



= 5.6 X 10^''(Ma/10-4Mo)" -^^ G 
MQyr"^)^'^-^] respectively, where 
Ma is the total accreted mass and Ma is the accretion rate. It is shown that the belt prevents 
north-south heat transport by conduction, convection, radiation, and ageostrophic shear. This may 
explain why millisecond oscillations observed in the tails of thermonuclear (type I) X-ray bursts in 
low-mass X-ray binaries are highly sinusoidal: the thermonuclear flame is sequestered in the magnetic 
hemisphere which ignites first. The model is also consistent with the occasional occurrence of closely 
spaced pairs of bursts. Time-dependent, ideal-magnetohydrodynamic simulations confirm that the 
equatorial belt is not disrupted by Parker and interchange instabilities. 

Subject headings: accretion, accretion disks — stars: magnetic fields — stars: neutron — stars: 
rotation — X-rays: bursts 



1. INTRODUCTION 

Thermonuclear (type I) X-ray bursts are observed from 
70 of the 160 lo w-mass X-ray binaries (LMX Bs) dis- 
covered to date l|Strohmaver fc BildstenI 120031) . They 
recur every few hours to days, when the accreted sur- 
face layer of the neutron star ignites and is incinerated 
by hydrogen and helium burning. Brightness oscilla- 
tions with millisecond periods are obs erved during ther- 
monuclear X-ray burs ts in 13 LMXBs l|Mimo et alJl2002t 
iPiro fc Bildstenl2004f) . They arise during burst onset be- 
cause the stellar photosphere is temporarily patchy while 
the thermonuclear flame spreads from its ignition point 
to cover the star. 

A surprising property of the burst oscillations is that 
they often persist into the tails of bursts^ after th e flame 
is exp ected to have engulfed the star ifStrohmavcr et aP 
Il997jl . Equ ally surprising is h ow sinusoidal the oscil- 
lations are. iMuno et al.l l)2002|) analyzed the harmonic 
content of 59 oscillations from six sources, observed with 
the Rossi X-Ray Timing Explorer (RXTE), and found 
that the Fourier amplitudes of integer and half-integer 
harmonics are less than 5 and 10 per cent of the maxi- 
mum signal respectively. These data imply that if there 
is one hot spot on the surface, it must lie near the rota- 
tional pole or cover an entire hemisphere, whereas if there 
are two, antipod al hot spots, they must lie near the ro- 
tational equator fMun o et alJl200^ . Most theories that 
explain why the oscillations persist into the tails of the 
bursts, e.g. uneven heating /coo ling during photospheric 
uplift ijStrohmaver et ahl Il997|) , or cyclones driven b y 
zonal shear in a geostrophic flow l|SDitkoyskv et all2002D . 
are hard pressed to account for the pattern of hot spots 
implied by the Fourier data. Global r-modes in the neu- 
tron star ocean can divide the photosphere into symmet- 

Electronic address: |dpayne@phy sics . unimelb . edu. au] 



ric halves l)Hev]ll20d^ . but the physical mechanism con- 
verting r-mode density perturbations to brightness oscil- 
lations is unclear (Muno et al. 2002). Gumming (200^ 
showed that differential rotation between the pole and 
equator of > 2 per cent can also excite unstable modes. 

In this paper, we show that an equatorial belt of intense 
magnetic field, compressed by accreted material spread- 
ing away from the magnetic poles, can impede thermal 
transport between the hemispheres of the star. In §2, 
we review the physics of magnetic burial and compute 
the maximum field strength and width of the equato- 
rial magnetic barrier. In §3, we estimate the efficiency 
of thermal transport across the barrier by conduction, 
convection, radiation, and ageostrophic shear, and in- 
vestigate whether cyclonic flows can disrupt the barrier. 
The implications for the harmonic content of burst oscil- 
lations are explored in §4. 

2. MAGNETIC BURIAL 

2.1. Grad-Shafranov equilibria 

In the process of magnetic burial, material accreting 
onto a neutron star accumulates in a column at the mag- 
netic polar cap, until the hydrostatic pressure at the base 
of the column overcomes the magnetic tension and mat- 
ter spreads equatorward, dragging along frozen-in polar 
magnetic field lines to form an equato rial magnetic belt 
or 'tutu'. FigureCJa), reproduced from lPavne fc MelatosI 
0004) , illustrates the equilibrium configuration obtained 
for A/a = 10~^Mo, where Ma is the total accreted mass. 
The polar mountain of accreted material (dashed con- 
tours) and t he pinched, flaring, equatorial magnetic belt 
are e vident ijMelatos fc Phinnevll200lt iPavne fc MelatosI 
|2004|) . The equatorial magnetic field strength increases 
in inverse proportion to the surface area of the equato- 
rial belt, by flux conservation. The Lorentz force per unit 
volume exerted by the compressed equatorial field [Fig- 
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ure^b)] balances the thermal pressure gradient [Figure 
nic)], and gravity, preventing the accreted matter from 
spreading to the equator. 

In the steady state, the equations of ideal magnetohy- 
drodynamics (MHD) reduce to the force balance equa- 
tion (CGS units) 



Vp + pS/(t> - (47r) (V X B) X B = 0, 



(1) 



where B, p, p = c^p, and 



GAUr/Rl denote the 



magnetic field, fluid density, pressure, and gravitational 
potential respectively, Cs is the isothermal sound speed, 
is the mass of the star, and i?* is the stellar radius. In 
spherical polar coordinates (r, 6, 0), for an axisymmetric 
field B — \/'i(j{r, 9) I {r sin 6) x e^, equation ^ reduces to 
the Grad-Shafranov equation 



A^^ = F'{'^) exp[-(0-0o)A 



where 
A2 = 



1 



fioT^ sin 9 



sin 9 d 



1 d 



ae \ sin 9 89 



(2) 



(3) 



is the Grad-Shafranov operator, F(ip) is an arbitrary 
function of the ma gnetic flux and we set 0o — 4'{R*)- 
In this paper, as in iPavne &: Melato^ pOO^ . we fix F(V') 
uniquely by connecting the initial dipolar state (with 
Ma = 0) and final distorted state [e.g. Figure ^^a)] via 
the integral form of the flux freezing condition of ideal 
MHD, viz. 

dM „ f dsn 

Here, C is any magnetic field line, and the mass-flux 
distribution dM/dip is chosen to capture the magneto- 
spheric geometry, e.g. the accretion stream is tunneled 
magnetically onto the pole, with dM/dip oc exp{—ip/ipi^), 
where tp;^ is the polar flux. We also assume north-south 
symmetry and adopt the boundary conditions "0 = dipole 
at r = i?* (line tying), ip — a,t 9 — 0, and dip/dr — 
at large r. The line- tying approximation artificially pre- 
vents accreted material from sinking, so the computed 
p is an upper limit. Equations (2) and (3) are solved 
numerically using an iterative relaxation scheme and an- 
alytically by Green functions, producing equilibria like 
Figure 

2.2. Equatorial magnetic belt 

The compressed magnetic field in the equatorial belt 
emerges approximately perpendicular to the stellar sur- 
face, with opposite sign in the two hemispheres. Near 
the equator, B falls off roughly as exp[— 0.7(7r/2 — 
6')/A6'] exp[-(r - R^)/h], where h = c'^RJGM^ is the 
hydrostatic scale height, A9 is the belt thickness, and 
the factor 0.7 comes from empirically fitting to the nu- 
merical results. The maximum magnetic field strength 
in the belt, -Bmax, computed numerically as a function of 
Ma, is plotted in Figure |21 It is fitted by 

f 2.0 X 10i6(Ma/Mc)°-9i±'' °6G Ma < 0.4Mc 

Bmax = < 6.3 X 10l5(Ma/Mc)°-32±O OlG Ma > 0.4Me 



with Mr = GM^BlRli 



IQ-^Mo 



(5) 

where Bq is the 



polar magnetic field strength prior to accretion. The 
scaling (0 agrees with analytic theory for small M^. 



(jPavne fc Melatosll2004|) . Also plotted in Figure His the 
half- width half-maximum thickness of the belt, fitted by 
M = 3°(Ma/Mc)-i-5±° °3. We find that M decreases 
as -Bmax increases, as expected from magnetic flux con- 
servation, but not exactly as A9 oc i?~axj because Br is 
underestimated numerically at the equator by ~ 10 per 
cent (flux loss due to finite grid resolution). Note that 
equation Q and the above scalings of A6' versus Ma do 
not include Ohmic diffusion, which is discussed further 
in §2.3 and g| 

Grad-Shafranov equilibria are difficult to compute di- 
rectly from (O and Q for Afa ^ 1.6Mc, because the 
magnetic topology changes abruptly: magnetic bubbles 
form which are disconnected from the surface, hinting 
at time-dependent processes which the steady-state the- 
ory cannot describe. In addition, the iterative relax- 
ation scheme struggles to handle steep field gradients. 
Therefore the results in Figure |2] for Ma > 1.6Mc are 
computed by direct numerical simulation using ZEUS- 
3D, a multipurpose, time-dependent, ideal-MHD code 
for astrophysical fluid dynamics which uses staggered- 
mesh finite differencing and operat or splitting in three 
dimensions ijStone fc N ormanll 19921) . We load the Grad- 
Shafranov equilibrium for Ma = 1.6Mc into ZEUS-3D; 
double Afa quasistatically over 250 Alfven times with 
inflow (and hence radial B) boundary conditions; stop 
the inflow and allow B to relax to a dipole at large r\ 
then iterate to reach Afa ~ 1O~^M0 l|Pavne fc MelatosI 
120061) . An isothermal equation of state is chosen and 
self-gravity is switched off. The experiment is performed 
for h/R^, — 2 X 10^^ (for computational efficiency) and 
scaled to neutron star conditions (/i/i?* = 5 x 10~^) ac- 
cording to Bmax oc this scaling is verified numeri- 
cally for the range of Ma in Figure El 

2.3. Stability and Ohmic relaxation 

Distorted ideal-MHD equilibria are often disrupted by 
Parker and interchange instabilities. Remarkably, how- 
ever, this is not true for the equilibrium in Figure ^ 
Figure 13 shows the results of an experiment in which 
the equilibrium is loaded into ZEUS-3D, perturbed, an d 
evolved for 500 Alfven times ijPavne fc MelatosI l2006|) . 
[The Alfven time is defined as h/v\, where va is the 
Alfven speed averaged over the grid.] It is marginally sta- 
ble; i?max oscillates via magnetosonic (phase speed « Cg) 
and Alfven (phase speed ~ ua) modes which are damped 
by numerical dissipation. The configuration already rep- 
resents the end point (reached quasistatically) of the non- 
linear Parker instability of an initially dipolar field over- 
laid with accreted material. It is not interchange unsta- 
ble because line tying prevents fiux tubes from squeezing 
past each other. That said, we caution that we have 
not yet investigated the full gamut of three-dimensional 
MHD instabilities; cross-field mass transport cannot be 
categorically ruled out. 

Another pathway to cross-field mass transport is if the 
magnetic belt relaxes by Ohmic diffusion, either dur- 
ing quiescence or during a burst itself. During quies- 
cence, diffusion occurs most rapidly at the base of the 
accreted layer. In the relaxation time approximation, 
with Coulomb logarithms set to 10, the electron-phonon 
conductivity for a hydrogen-helium mixture is cr = 6.3 x 
10^''(/9/10^^g cm~^) s~^, the density at the base of the 



reted layer is p ^ 6.2 x lQ'^°{M^/lQ-^MQf/'^ 



g cm 
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ijBrown &: BildstenI 11998(1 . and hence the Ohmic diffu- 
sion time-scale across the equatorial belt (at the base 
of the accreted layer) is = AiraRlAO'^ /c^ = 2.6 x 

10^(Ma/10-5MQ)-9/4(y/208K)-iyr, much longer than 
the burst time-scale, where we rewrite A9 in terms of 
Afa- During a burst, diffusion occurs most rapidly in 
the burning layer. The temperature of the burning layer 
rises isobarically until the radiation pressure dominates 
the hydrostatic pres sure, reachin g T ^ 10* K for a typi- 
cal ignition column l|Brownll2004|) . The elevated temper- 
ature and reduced density lower a and accelerate Ohmic 
diffusion. In the burning layer [p « 10^ g cm~'^), we 
find td — 1.5 X 10^ yr, still much longer than the burst 
time-scale. 

Ohmic diffusion will be modelled self-consistently in a 
future paper. 

3. THERMAL TRANSPORT ACROSS THE EQUATOR 

Does the equatorial magnetic belt impede thermal 
transport enough to stop the thermonuclear flame in 
a type I X-ray burst from spreading from one hemi- 
sphere to another? A burst is initiated locally by 
a thi n shell thermal instability l|Schwarzschild fc HarmI 
Il965j) . As the nuclear burning time-scale is much shorter 
than the time to accrete the minimum column for igni- 
tion, the accreted layer ignites at a single point, most 
likely at the equator where gravity is rotationally re- 
duced, and the thermonuclear flame spreads away ei- 
ther as a def lagration front ( Fryxell fc Wo oslev 198^ 
Bildstenlll99a). bv detonation ifFrvxell fc Woo slev 1982; 
ZinEfale et al.l200lH . or as a cyclone driven by zonal shear 
ijSpitkovskv et alJ20 02;l. Detonation, which occurs when 
the nuclear burning time-scale is less than the vertical 
sound crossing time-scale, requires a thick 100 m) col- 
umn of fuel and hence a low accretion rate (< 10~"'^^'^Mq 
yr""'^). Deflagration occurs most commonly, with the 
front propa gating at a speed set by the heat flux and 
convection l|Frvxell fc Wooslev"1982^. Note that tem- 
perature fluctuations leading to single-point ignition are 
difficult to create because the sound crossing time around 
the star is very short, but Coriolis forces can assist by 
balancing sidew ays pressure gradients w hen the star is 
rapidly rotating ISpitkovskv 

We estimate below how the equatorial magnetic belt 
modifies heat transport in these scenarios. 

3.1. Conduction 

The thermal conductivity k of degenerate electrons at 
the atmosphere-ocean boundary of a magnetized neu- 
tro n star, where H /He burning occurs, was calculated 
by iPotekhinI l|1999() . We extract approximate values 
for the conductivity perpendicular to the magnetic field 
from Figures 5 and 6 of that paper, obtaining k±_ « 
10^-5(/9/10^gcm-3) (T/10*K)2(B/1015G)-2 in units of 
erg s~^cm~^K~^ in the regime 10'' g cm^'^ ^ P ^ 10^ 
g cm-3, 10^2 G < B < 10^5 G, and lO*' K < T < 10* 
K. Note that k± at the equator is reduced 10^ times rel- 
ative to its value before accretion commences (oc 
and, for reference, one has k±/k\\ ~ 10~^ aX T — 10* K, 
B = 10^2 G llPotekhinlliggl . We estimate the conduc- 
tion time-scale across the magnetic barrier from icond — 
RlA0^Cp/K±, where C = 10''(T/10*K) erg g-^ K"^ is 
the specific heat capacity of a degenerate Fermi gas 



(|Brown fc Bildstenlll99^ . finding 

tcond= 27(Ma/Mc)-3(T/10*K)-i(B/10i5G)2yr, (6) 

which is safely longer than the duration of the burst. 
Note that this estimate does not include corrections due 
to electron-electron scattering and impurity scattering. 
Furthermore, our fit smoothes over the wiggles in the 
solid curves of Figures 5 and 6 in IPotekhinI l|1999(l . which 
arise from quantization into Landau orbitals. 

3.2. Radiation 

In a strong magnetic field, photons polarized perpen- 
dicular to B dominate radiative transport. In the burn- 
ing (helium) layer, where Thomson scattering dominates, 
the Rosseland mean opacity _L B is given by = 
1.3 X 10-S(r/10*K)2(S/10i'"^G)-2g-i cm2, decreasing in 
inverse prop ortion to the square of the cyclotr on fre- 
quency {lOc ) (|.Toss fc LilllflMlFrvxe]] fc Wooslevl h 9821. 
Hence the optical depth of the barrier is tr = KRpi?* A0, 
i.e. 

TR = 6.8 X 10^ (Afa/Mc)~3/2(y/;L0*K)2(S/10i'"^G)-^ 

The optical depth Q is a lower limit obtained by assum- 
ing that radiation is transported from one hemisphere to 
the other near the surface, at the depth of the burning 
layer {p ^ lO^gcm"^) rather than at the base of the ac- 
creted column {p ^ lO^^gcm"'^). Hence the peak flux 
penetrating the barrier is Fburst — -^burste~^^/47ri?^ — 
8 X lO^'^e"^'^ erg s'^ cm'^ (for Lburst = lO^^ergs'^), 
which is insufficient to heat the other (quiescent) hemi- 
sphere above its ignition temperature 10* K) (as- 
suming thermal equilibrium and applying the Stefan- 
Boltzniann law). Interestingly, once -Bmax exceeds 10^^ 
G, we find tr < 1 and the magnetic belt becomes op- 
tically thin. Vertical heat propagation is not considered 
here. 

Equation {71) can be generalized in several ways. Mode 
coupling (Miller 199jl) can cause a fraction ~ 0.2uj/uJc of 
the perpendicular mode photons to convert into parallel 
mode photons, so that the opacity scales oc B~^, increas- 
ing the field strength required for the belt to become opti- 
cally thin. Vacuum polarization ('Ozer'2003), along with 
proton-cyclotron resonance, gives sharp spikes in the fre- 
quency response of the atmospheric opacity. However, 
the resonant densities for vacuum polarization (p ~ 10"'^ 
g cm~^) are achieved at shallow depths, well above the 
burning layer. 

3.3. Convection 

iBildsterJ 11199511 estimated the convective speed in 
the burning layer in terms of the mixing length, 
Im, and thermal time-scale, tth, finding Vc ~ 
Csilm/h)^^^{csRl/GMJthV^^ « lO^cms-i. This is 
consistent with the upper bound 10 ''cms""' obtained if 
^burst is transported entirely by the mechanical flux pv^ . 
It is also co nsistent with the obser ved burst rise time 
R*/vc ^ 1 s (|Spitkovskv et alJl2002|) . Convection is sta- 
bilized magnetically if the magnetic tension exceeds the 
ram pressure pv^, which occurs for Smax ^ (STrpf^)'^^ ~ 
5 X 10^ G. This condition is met comfortably in the equa- 
torial belt^^lie_jM£neticfi can also quench convec- 
tion l|Gough fc Tavleil[T96l . 
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3.4. Ageostrophic shear flow 

ISpitkovskv et al.1 1)20021) suggested that inhomogeneous 
cooling drives zonal currents which are unstable to the 
formation of cyclones, as in planetary atmospheres. This 
may explain why the coherent oscillations observed in the 
tails of some type I X-ray bursts persist for many rise 
times. The drift in oscillation frequency Hz) during 
the burst is attributed to the Coriolis drift of the cyclone 
in the frame of the star, although theory predicts larger 
frequency drifts 10 Hz) than those observed. The os- 
cillation amplitude in the burst tails (~ 10%) is governed 
by processes other than magnetic burial, e.g. small-scale 
magnetic fields generated by an MHD dynamo in the 
burning front, which are confined to the ashes after the 
burst, w hile the freshly accreted matter remains unmag- 
netized l)Spitkovskv et alJl2002|) . 

Ageostrophic shear flow moves hot mate- 
rial ahead of the burning front and draws 
fresh fuel into it at the flame speed inflame ~ 
2x10^ cm s-i(/ihot/103cm)(//0.32kHz)-i(W/0.1s)-i, 
where /ihot is the scale height of the incinerated ocean, 
/ — 2ilcos9 is the local Coriolis parameter, is 
the angular frequency of the star, and t^uc is the 
nuclear burning time-scale. Magnetic tension stabilizes 
ageostrophic shear for B > {STrpvl^^^J^/'^ w 4.5 x 10^" 
G <^ -Bniax- In other words, when the cyclonic flame 
runs into the equatorial magnetic belt, it is reflected; 
conversely, ageostrophic shear cannot disrupt the belt. 
Latitudininal shea r instabilit i es, a possible source of 
burst oscillations l) Cummin i |2005), may disrupt the 
magnetic belt and provide one motivation for extending 
our burial model to three dimensions in the future. 

4. DISCUSSION 

Polar magnetic burial creates an intense, equatorial 
belt of magnetic field which can thermally isolate the 
magnetic hemispheres of an accreting neutron star. The 
maximum magnetic field strength in the belt, -Bmax — 
5.6 X 10^"^{MjlO-'^MQf-^^ G, is sufficient to prevent 
heat transport by conduction, radiation, convection, and 
ageostrophic shear. The conduction time-scale tcond 
27 yr exceeds the cooling time of the incinerated material; 
the magnetic belt is opaque (optical depth 7 x 10^); 
convection is stabilized magnetically (pu^ Bf^^^/8n); 
and so is ageostrophic shear {pvg^^^^ ^ B^^^^/8n). How- 
ever, the conclusion that the hemispheres are thermally 
isolated by the magnetic belt is less secure at large ac- 
creted masses (Ma ~ O.IM©), where the Grad-Shafranov 
and ZEUS-3D calculations in section 2 are hampered by 
numerical difhculties, and Ohmic diffusion (which we do 
not incorporate self-consistently) becomes important. 

Thermal isolation of the magnetic hemispheres is con- 
sistent with the highly sinusoidal light curves of mil- 
lisecond oscillations in thermonuclear bursts in LMXBs 
(^Imio et al. 2002). During the rise of the burst, oscilla- 
tions are caused by the spre ading of a hot spot, pr oba- 
bly in the form of a cyclone l)Spitkovskv et alJl2002D . In 
the tail of t he burst, the hot spot is sequestered in one 
hemisphere ()Muno et alJl200^ . and misalignment of the 
magnetic and spin axes guarantees that we observe per- 
sistent oscillations at the spin frequency. Furthermore, it 
is observed that bursts occasionally occur in quick suc- 
cession, separated by a 5 — 10 minute interval, compared 



to an interval of several hours between typical bursts 
l)Lewin et al.lll993l) . Such burst pairs are consistent with 
the model of magnetic burial: if a burst ignites one hemi- 
sphere, the other remains dormant and can, in principle, 
ignite shortly thereafter when more material accretes. If 
this explanation is valid, one would expect that there are 
no triple bursts, that the fluences emitted in each half 
of a burst pair are roughly equal, and that the number 
of burst pairs relative to normal bursts increases as the 
ignition probability per unit time increases (relative to 
Ma).^ A more detailed calculation is required to resolve 
whether the 10-min intervals are caused by delayed heat 
propagation through the magnetic belt. 

Muno et al. (2002) showed that the bright burning re- 
gion must cover 80° — 110° in latitude in order to match 
the upper limit on the observed ratio of harmonic to 
fundamental amplitudes, which implies A9 < 10°. If 
the scaling of A9 versus Ma in i)2.2l is taken at face 
value, this requires Afa > 0.5Mc. Furthermore, for 
the belt to be optically thick, we need Ma < 2Mc 
f H3.2|l . Apparently then, only a narrow range of accreted 
masses (0.5 < M^/Mc < 2) can account for the obser- 
vations. However, the above A6 scaling, which depends 
steeply on Ma, does not tell the whole story, because 
Ohmic diffusion is not incorporated self-consistently. 
Equatorward hydromagnetic spreading is arrested when 
the accreti on time-scale exceeds t he Ohmic diffusion 
time-scale l)Brown fc BildstenI 11998)) . softening the de- 
pendence o f i?ma,x and AO on M„ in the manner de- 
scribed by iMelatos fc Pavnel ((2005). Ohmic diffusion 
arrests magnetic compression for A/a > Afd, where 
Mh = 3.4xlO-'^(r /10^K)-'^-^(A/a/10-^Mr^vr-M°-'^^M^ 
l)Melatos fc Pavn'3 120051) is the accreted mass at which 
the accretion time-scale exceeds the Ohmic diffu- 
sion time-scale. Including this effect, A6 is mod- 
ified to the maximum of 3°(Ma/10-''MQ)-i-^ and 
3°(Afa/lO-''Mo)O-5(A/a/lO-*M0yr-i)-O-5], and the belt 
remains optically thick even for Afa ^ 2Mc f ti3.2|) . 

Only 70 out of 160 LMXBs undergo bursts, and only 
13 exhibit millisecond oscillations. Does the equatorial 
belt model respect these statistics? If accretion occurs 
at the Eddington rate, A/a ~ 10~®Moyr~^, it takes less 
than 10'* yr to achieve Ma > Mc and screen the po- 
lar magnetic field, allowing bursts to ignite and creat- 
ing a thermally insulating equatorial belt. However, for 
Afa ^ 2 X lO^Afc, the belt becomes optically thin to X- 
rays and tcond drops below the duration of the burst, 
potentially allowing the flame to engulf the entire star 
in LMXBs older than a few times 10^ yr. If accretion 
occurs at Aifa ^ 10~^^Mp,yi~^ , as in a ccreting millisec- 
ond pulsars ()Chakrabartv et all I2003D . it takes 10'' yr 
to achieve Ma > Mc, screen the polar magnetic field, 
and allow bursts to ignite. Therefore, we do not ex- 
pect bursts from all accreting millisecond pulsars, but 
when bursts do occur, we expect to detect oscillations at 
some level because we have tr ^ 1 and tcond ^ (burst 
time-scale) for Afa < 1O~^M0. Such is the case, within 
current observational sensitivity, for SAX J1808. 4-3658 

^ If the ignition probability is higher (given a specific Ma), it 
is more Ukely that the second hemisphere will ignite shortly after 
the first (before the first has a chance to refill and ignite again). A 
faster burst rate may be associated with more pairs, but gaps in 
the satellite data make it difficult to be sure. 
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and XTE J1814-338. Sources like XTE J1814-338, whose 
oscillations contain harmonics exceeding 25 per cent of 
the peak amplitude ijStrohmaver et ani2003j) . may have 
semi-transparent equatorial belts; perhaps this millisec- 
ond pulsar has experienced interruptions in its accre- 
tion history resulting in Ma < Mc- Note that, for 
-^^a ^ 10^^°Moyr~^, Ohmi c diffusion may magn etize 
the freshly accreted material i Cumming et al.ll2001fl . in- 
creasing the polar magnetic field strength and suppress- 
ing bursts. Also, ignition is more likely near the equator, 
where gravity is rotationally reduced, but this is where 
the magnetic field is strongest. The i ssues of conser- 
vative mass-transfer (jTauris et al.ll20 0'?l and differential 
rotation leading to shear instabilities (Gumming 2005) 
are not considered here. 

On the face of it, high-mass X-ray binaries (HMXBs) 
accrete many times Mc, yet a strong magnetic field sur- 
vives. However, one needs to be cautious. While the 
accretion rate in HMXBs can reach f O~^'^-lO~^M0yr~^ 
during either atmospheric Roche lobe overflow (close bi- 
naries with orbital period < 5 days) or stellar wind accre- 
tion (wider binaries), the typical life time of HMXBs a s 
strong X-ray sources is fO -fO^ yr ljUrpin et alJlf99^ . 
Sub-Eddington accretion rates (lO~^^-lO~^°M0yr~^) 
characterize the rest of the companion's main-sequence 
evolution for 10^-10^ yr. These scenarios yield roughly 
Afa ~ Mc. Mass transfer can also be nonconservative, 
further reducing Ma, e.g., in intermediate-mass X-ray bi- 
naries (jTauris e t al. 200(|). Of course, despite these cau- 
tionary remarks, it may well be that some HMXBs do 
accrete many times Mc. If so, then either HMXBs are a 
counterexample to the simple magnetic burial model we 
have calculated, or else cross-field transport by Ohmic 
diffusion (which we do not incorporate self-consistently) 
becomes important. 
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A recent model for doubly peaked bursts which 
are too weak to cause p hotos pheric expansion 
(jBhattacharvva fc Strohmaved I2006D proposes that 
a burning front forms quickly after ignition at or near a 
pole and propagates quickly towards the equator where 
it stalls. After a delay, the burning front speeds up again 
into the opposite hemisphere. While the time between 
burst peaks (~ 10 seconds) is less than that estimated 
for burst propagation across the magnetic belt, a weak 
magnetic belt may explain why the burning front stalls 
at the equator. 

Numerical difficulties associated with steep field gradi- 
ents prevent us from verifying whether our model scales 
up to Ma ~ O.1M0. Howeve r, if we continue to respect 
mass-flux conservation as in ii2.1l it is probable that ac- 
creting extra mass does not eliminate the magnetic bar- 
rier. Even if, for example, there are instabilities which 
disrupt the equatorial magnetic belt above ~ 1O~^M0, 
the belt does not stay disrupted; magnetic burial ensures 
that it r eforms as soon as a further ~ 1O~^M0 is accreted 
ijPavne fc Mel atos 2004). Given that O.IM0 is typically 
accreted in LMXBs, the chances of catching the belt in 
its disrupted state are slim. 
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Fig. 1. — {top) Equilibrium magnetic field 
i/'a = O.lip*. Coordinates measure altitude. 
0.8, 0.6, 0.4, 0.2, 10-2, 10-3, IQ-", IQ-'^, IQ-'', 10 



lines (solid curves) and density contours (dashed curves) for M; 
Density contours are drawn for Typmax (Pmax — 2.52 X lO^"* 
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gem— ''), with ri = 
MelatosI 120041) . 1 (mid- 



dle) Contours of Lorentz force per unit volume for the same r] values, (bottom) Contours of pressure gradients for the same rj values. Note 
that colatitude = 0° at the pole. 
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Fig. 2. Maximum magnetic field in the equatorial belt, _Bmax, computed numerically for = 2 X 10 ^ (crosses) and h / Rt = 5 X 10 ^ 
(triangles) and scaled using -Dmax oc h ^, plotted together with the half-width half-maximum thickness, A6 (squares), as a function of 
accreted mass, Ma. 
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Fig. 3. — Maximum magnetic field strength, -Bmax, as a function of time (in units of the Alfven time) for M^/Mc = 
0.16, 0.32, 0.64, 1.12, 2.4, 4.0 (bottom to top) when the equihbrium in Figure is loaded into ZEUS-3D and perturbed slightly. The 
configuration is stable. 



